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摘 要 : 在 太阳 观测 研究 中 ， 高 分 辨 图 像 与 全 日 面 像 的 配 准 是 一 项 非常 有 意义 的 工作 ， 
但 由 于 他 们 之 间 存 在 着 旋转 、 缩 放 和 平移 ， 因 此 很 难 高 精度 地 进行 匹配 。 提 出 一 种 结合 局 部 
统计 信息 和 控制 点 匹配 的 图 像 配 准 方法 , 核心 思想 为 将 视 场 等 间隔 划分 为 大 量 重 登 的 局 部 区 
域 , 通过 相关 匹配 在 全 日 面 像 上 寻找 对 应 的 局 部 区 域 , 然后 计算 每 一 对 局 3 部 区 域 间 的 亚 像 元 
偏 移 ， 根 据 偏 移 量 确定 每 一 对 特征 点 的 坐标 位 置 ， 据 此 作为 点 匹配 中 的 特征 控制 点 ; 最 后 根 
据 控 制 点 建立 仿 射 变 换 的 转换 方程 , 采用 最 小 二 来 求解 整个 视 场 的 转换 参数 ， 生出 的 人 
数 重新 对 图 像 进行 选 代 ， 得 到 收敛 后 的 结果 并 进行 配 准 。 通 过 对 高 分 辨 观测 图 像 和 全 
SDO/HMI 连续 谱 图 像 进行 配 准 ， 拟 合 结果 的 偏差 在 0.25” 以内。 
关键 词 : 太阳 图 像 配 准 ; ”高 分 辨 图 像 ; 全 日 面 图 像 ; 
中 图 分 类 号 : P182. 2 ”文献 标识 码 : A 文章 编号 : 1672-7673 


在 研究 太阳 图 像 的 观测 数据 时 , 经 常 需要 将 不 同 设备 的 观测 数据 进行 对 比 ， 从 而 发 挥 它 
们 各 自在 不 同 波段 、 不 同 视 场 、 不 同 分 辨 率 上 的 观测 优势 。 太 阳光 学 成 像 观 测 包 括 全 日 面 和 
高 分 辨 观测 两 类 。 对 于 全 日 面 图 像 ， 由 于 太阳 几何 尺寸 和 运动 方式 已 知 ， 如 太阳 的 半径 、P0 
角 、 自 转 等 ， 因 此 各 种 观测 设备 均 把 自己 的 数据 在 预 处 理 后 归 算 到 标准 的 太阳 坐标 系 。 然 而 
XIF fae mers 由 于 观测 设备 的 指向 误差 、 视 场 旋转 以 及 焦距 的 变化 ,使 得 直接 依赖 仪 
器 参数 确定 准确 的 太阳 视 场 坐标 是 不 精确 的 , 因此 只 能 借助 太阳 图 像 上 相似 的 太阳 结构 (如 
光 球 层 对 应 光 球 层 ， 色 球 层 对 应 色 球 层 ) 与 全 日 面 像 进行 转换 ， 即 视 场 匹配 。 然 而 全 日 面 观 
测 有 大 于 31” 的 视 场 , 但 最 高 的 空间 像 元 分 辨 仅 为 0.5”( 空 间 太 阳 动 力 天 文 台 的 日 球 层 磁 
场 观测 仪 ，SDO/HMI)。 而 地 面 高 分 辨 观测 ， 如 中 国 的 lm 新 真空 太阳 望远镜 (New Vaccum 
Sloar Telescope,NVST) 和 美国 大 熊 湖 天 文 台 的 新 太阳 望远镜 (New Solar Telescope,NST), 视 
场 仅仅 为 1” 到 2”， 但 像 元 分 辨 率 达 到 0.05” 以 上 ， 这 就 意味 着 需要 将 两 个 空间 分 辩 率 相 
差 十 倍 的 图 像 进行 配 准 。 即 便 不 考虑 其 他 畸变 ， 局 部 像 的 视 场 旋转 CRD. EARN CS) 和 中 
心 位 置 坐标 “7x，7y) 是 必须 得 到 的 4 个 参数 。 
图 像 配 准 是 指 在 不 同时 间 、 不 同 条 件 下 得 到 的 并 且 位 于 不 同 坐 标 系 下 的 同一 场景 的 两 幅 
或 者 多 幅 图 像 进行 匹配 的 过 程 "。 由 于 存在 各 种 不 同 的 情况 ,同一 场景 的 多 幅 图 像 会 在 分 辨 
率 、 位 置 、 尺 度 、 成 像 模 式 等 方面 存在 很 多 差异 ， 图 像 配 准 是 找到 一 个 几何 变换 的 最 优 解 ， 
使 得 图 像 上 相对 应 的 点 , 或 者 至 少 是 有 意义 的 点 达到 [匹配 状态 。 目 前 在 图 像 配 准 领域 , 常见 
的 配 准 方 法 有 基于 特征 的 点 匹配 方法 , 即 根 据 图 像 中 的 特征 点 , 或 者 显著 特征 构造 描述 子 确 
定 的 匹配 点 ， 然 后 构造 转换 方程 进行 配 准 。 理 论 上 这 个 方法 如 果 有 足够 多 和 足够 精确 的 控 带 
点 是 可 以 很 好 地 解决 两 个 图 像 的 几何 变换 ， 如 旋转 、 缩 放 、 平移 等 问题 。 但 关键 问题 是 如 何 
确定 这 些 点 ， 这 方面 在 图 像 处 理 界 有 很 多 研究 ， 如 Harris 角 点 及 其 衍生 的 Harris-Laplace 方 
法 ”“， 尺 度 不 变 特征 (Scale-invariant feature transform,SIFT) “及 其 生生 的 PCA-SIFT 算法 
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(Principal Com-ponent Analysis SIFT)"、 局 部 特征 描述 子 ”、 高 速 尺度 不 变 特征 算法 (Very Fast 
SIFT, VF-SIFT) “， 加 速 鲁 棒 性 特征 (Speeded Up Robust Features,SURF)““ 等。 另外 一 
种 是 基于 区 域 信息 的 统计 配 准 方法 , 主要 通过 图 像 区 域 的 统计 信息 度量 图 像 的 相似 程度 ， 如 
互信 息 ”， 互 相关 函数 "”"， 模 板 匹 配 “"， 不 变 矩 ， 相 关 匹 配 *”" 等 。 但 一 般 来 讲 这 种 方法 
主要 解决 图 像 间 的 平移 问题 , 对 于 具有 旋转 和 比例 缩放 的 图 像 ， 又 引申 出 在 频率 域 的 相位 相 
eR, Fouier-Mellin "等 。 
然而 ， 由 于 太阳 观测 图 像 和 普通 自然 图 像 有 很 大 的 区 别 ， 其 中 的 特征 边界 不 清晰 ， 结 构 
相对 模糊 且 存 在 很 多 噪声 , 同时 动态 范围 大 , 运动 速度 较 快 ,并且 存 在 很 多 不 同 的 活动 现象 ， 
大 部 分 运动 现象 为 非 刚 性 运动 ; 同时 太阳 图 像 的 图 幅 一 般 较 大 ,有 效 目标 小 且 缺 少 有 效 的 纹 
H, 不 存在 明显 的 梯度 变化 等 特征 , 并 且 由 于 拍摄 太阳 图 像 的 设备 以 及 拍摄 过 程 中 的 不 确定 
PE, 采集 的 太阳 图 像 还 存在 旋转 、 缩 放 及 平移 等 情况 , 不 同 图 像 的 空间 分 辩 率 经 常 不 在 同一 
尺度 ， 这 就 为 确定 特征 对 应 点 带 来 非常 大 的 困难 。 长 期 以 来 , 在 同一 观测 设备 的 太阳 时 间 序 
列 图 像 的 视 场 对 齐 一 直 采 用 基于 互相 关 的 统计 配 准 ， 而 不 同 观 测 设备 之 间 的 图 像 配 准 ,还 停 
留 在 利用 一 些 显著 的 太阳 特征 ， 如 黑子 、 气 孔 、 大 米粒 构造 对 应 点 的 方式 。 然 而 这 些 特 征 点 
g 数量 少 ,结构 不 规则 ,难以 保证 位 置 测量 精度 , 同时 也 难以 使 用 计算 机 技术 进行 自动 化 测量 。 
~ 尤其 是 将 高 分 辩 局 部 像 对 应 到 全 日 面 像 的 时 候 ， 这 种 问题 更 为 突出 。 但 从 另 一 方面 , 通过 望 
© 远 镜 机 械 和 终端 光学 系统 ， 对 局 部 太阳 像 指 向 、 比 例 尺 和 位 置 有 初步 估计 ， 只 是 这 些 数值 不 
gs 很 精确 且 随 时 间 有 一 些 变化 。 因 此 实际 面临 的 是 如 何 提高 配 准 精度 ,进行 精细 图 像 匹 配 的 问 
© 基于 这 些 问 题 ， 本 文 提 出 了 一 种 将 统计 信息 和 相关 点 匹配 相 结合 的 太阳 图 像 配 准 算法 ， 
T 解决 将 不 同 分 辨 率 下 的 局 部 高 分 辨 观测 数据 和 全 日 面 数据 进行 高 精度 坐标 转换 的 问题 , 其 中 
包括 测量 图 像 指 向 角 、 像 元 比例 尺 和 图 像 中 心 在 全 日 面 像 上 的 位 置 。 这 一 方法 的 基本 思想 是 
在 初步 匹配 的 基础 上 , 首先 将 视 场 划分 成 大 量 重 登 的 局 部 区 域 代替 点 匹配 中 的 控制 点 ,通过 
基于 统计 的 相关 匹配 确定 对 应 区 域 间 的 亚 像 元 位 移 , 从 而 确定 特征 点 对 在 同一 坐标 系 中 的 位 
置 ， 然 后 用 最 小 二 乘 迭 代 求 解 整个 视 场 的 转换 参数 ， 并 逐步 迭代 提高 拟 合 的 精度 ， 成 功 将 
IM 太阳 望远镜 和 新 太阳 望远镜 的 Ti0 观测 数据 与 几乎 同时 刻 《〈 最 长 间隔 11 s) 的 SDO/HMI 
连续 谱 进 行 了 匹配 ， 拟 合 结果 的 偏差 可 达到 0. 25” 之 内 。 
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1 实验 数据 来 源 


本 文 实验 使 用 的 全 日 面 图 像 来 自 太阳 动力 学 观测 卫星 SDO(CSolar Dynamics Observatory, 
SDO) 上 日 球 层 磁 场 观 测 仪 HMI (Helioseismic and Magnetic Imager, HMI) 采集 的 连续 谱 全 
日 面 像 ， 而 高 分 辨 图 像 为 IM 太阳 望远镜 和 新 太阳 望远镜 的 TiO 光 球 像 。SDOVHMI 设备 
的 空间 分 辨 率 约 为 0. 5043 arcsec/pixel， 图 像 尺 寸 为 4096 像 元 *4096 像 元 ，1M 太阳 望 远 
镜 和 新 太阳 望远镜 图 像 的 像 元 比例 尺 为 0. 04 arcsec/pixel 和 0. 03 arcsec/pixel。 

图 1 Ca) 为 1M 太阳 望远镜 的 Ti0 观测 数据 ， 观 测 时 间 为 2013-7-15 07:29:09， 有 效 
区 域 尺寸 为 2304 像 元 *1920 像 元 ; 图 1 (b) 为 SDO/HMI 连续 谱 观 测 数据 ， 观 测 时 间 为 
2013-7-15 07:29:15, 红 框 标记 区 域 为 IM 太阳 望远镜 在 SDO/HMI 中 的 大 体位 置 。 本 文 主 
要 以 此 图 像 的 匹配 说 明 算 法 流程 。 
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图 1 观测 数据 


Fig. 1 Observed data 


2 配 准 的 方法 


2. 1 配 准 方法 的 流程 
配 准 方法 的 主要 流程 如 图 2。 
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图 2 算法 流程 
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Fig.2 Algorithm flow chart 
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计 ， 得 到 旋转 (Rotation). Ai 


图 2 可 以 描述 为 下 列 11 个 步骤 ; 
C1) 对 1M 太阳 望远镜 像 的 指向 、 比 例 尺 和 位 置 进行 初步 估 

放 (Scaling) KFE (Translation) 参数 (下 文 称 RST 参数 ); 
对 SDO/HMI 全 日 面 像 进行 旋转 ， 并 截取 与 IM AE 


望远镜 像 视 场 相 同 


(2) 根据 RST BR, 


的 SDO/HMI 子 图 太 ; 
(3) 根据 RST 参数 和 SDO/HMI 的 比例 尺 , 将 1M 太阳 望远镜 像 缩 小 至 和 SDO/HMI 同一 比例 尺 


下 ， 即 对 IM 太阳 望远镜 图 像 进行 降 采 样 ; 
(4) 对 降 采样 后 的 IM 太阳 望远镜 像 进行 高 斯 平滑 , 减 小 采样 过 程 中 的 频率 混 县 。 平滑 后 的 
H Tp RIN; 


IM 太阳 望远镜 降 采 样 像 有 
将 图 像 玉 等 间隔 划分 为 N 个 重 倒 的 子 块 ， 并 将 每 个 子 块 记 为 广 
得 大 对 应 的 每 一 个 到 


(5) 
(6) “对 每 一 个 子 块 产 和 SDO/HMI 子 图 到 利 用 归 一 化 互相 关 方法 求 得 
FH fis 

(7) “对 每 一 组 天 和 天， 测量 其 亚 像 元 偏 移 ， 并 将 子 块 中 心 设置 为 一 组 对 应 的 控制 点 ; 
(8) 筛选 每 组 控制 点 ， 将 控制 点 间距 离 大 于 一 定 像 元 的 点 剔除 ; 
7m (9) 根据 剩 下 的 控制 点 ， 建 立 仿 射 变换 的 转换 方程 ， 采 用 最 小 二 乘 求 解 RST 参数 ; 
~ (10) 计算 残 差 MRRKML, Hw EL <05^ <01}; { <03}; 
N {L>30, =min( ,_,,)} Et WRR: E, REPR (2)， 根 据 步骤 (9) 解 得 

行 迭代 处 理 ; 


| SDO/HMI 和 1M 太阳 望远镜 像 进 


的 RST 参数 重新 对 
(11) 最 终 得 到 迭代 完成 后 的 RST 参数 并 进行 配 准 
法 的 核心 步骤 ， 下 文 对 该 部 分 进行 更 为 


(6) 至 (9) 步 为 算 ; 
的 RST 参数 会 根据 上 一 次 迭代 的 结果 进行 下 一 


\ 体 描述 。 
Wik 


在 上 述 步 又 中 ， 

= mete, Fanta (2) 到 (3) 步 中 

q 代 ， 意 味 着 算法 会 重新 对 图 像 进 行 平 滑 。 
2. 2 分 块 确定 特征 控制 点 

EKR, 记 为 集合 {fi,i=1,2...,N}; 


设 冯 为 待 匹配 图 像 ,5 为 参考 图 像 。 将 Lr 划分 N 个 特 和 
六 的 尺寸 为 太 V ,每 个 区 域 有 的 尺寸 为 u v ， 可 知 M<UV<Y ， 则 参考 图 像 1; 和 图 像 fi 


的 互相 关 函 数 Corr(x,y) 定 义 为 ， 
SSe) Fale) F] 
Corr(x,y) = — t=0 一 一 一 一 (1) 
L (erst i] »> KTEJ 
A Reads ‘Ti (x+s,y+8)] ; A= ASe]. LA 


tH, (s:t) 0 fi ERBER fi = 


ZN 
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/和 fg Shi A ASS EX IRAP, RAL Ty A fi 覆盖 的 (5,1) 进行 


根据 上 述 公 式 ， 每 个 六 和 到 做 归 一 化 互相 关 ， 可 得 到 一 个 互相 关 曲 面 ， 而 峰值 的 坐标 


代表 了 fi 图像 中 心 在 1; 中 的 整数 像 元 精度 的 匹配 位 置 。 


在 1; 中 根据 上 述 位置 可 截取 和 友 同 样 大 小 的 特征 区 域 , 得 到 每 一 个 fi 对 应 的 大 。 由 它 


们 组 成 的 集合 记 为 {fx,i=1,2...,N}。 


对 于 两 组 对 应 的 特征 区 域 的 fi 和 fs， 使 用 集合 五 描述 六 和 Se 任意 两 个 特征 区 域 之 间 
的 关系 ， 即 


P=F(fi,fi)lsisN , (2) 


F, 表示 任意 两 个 对 应 的 局 部 特征 区 域 (f/,f: )。 


对 于 每 一 对 特征 区 域 ( 太 ,所 ) ,计算 标准 的 互相 关 函 数 。 假设 DD 为 特征 区 域 (fi, fi ) 间 


互相 关 函 数 分 布 的 峰值 局 部 区 域 ， 按 照 修正 矩 方 法 计算 重心 坐标 (m,n) , 


gz A) _ 2,0) 
= ee * 


(3) 


其 中 ，1/ 为 曲面 强度 值 ， 则 7'= 了 一 b, b=min(7),TeD。 i—i (fi, fi) 图 像 间 互 相关 函数 的 


局 部 重心 坐标 代表 每 一 对 特征 区 域 ( 太 ,fi ) 间 的 亚 像 元 偏 移 。 


根据 亚 像 元 偏 移 的 情况 ， 可 以 建立 每 一 对 ( 大, 大 ) 图 像 中 心 点 对 应 的 坐标 , 得 到 和 ,之 
间 的 NW 对 坐标 转换 控制 点 。 用 集合 P 描述 控制 点 对 之 间 的 关系 ， 即 


P=P( pp ps) 1Si<SN . (4) 


2. 3 最 小 二 乘 求解 转换 参数 
考虑 平移 、 旋 转 和 缩放 的 情况 ， 假 设 控制 点 对 (ph, Ps) ， 对 应 的 坐标 分 别 为 (X4,y4 ) 和 
(xs 加) ， 变 换 关 系 为 


x, =Ax+(1+m)x,cos@—-(1+m)y, sing 


(5) 
y, =Ay+(1+m)x, sinO+(1+m) y, cosO 


转换 为 矩阵 形式 ， 
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Ax 
x| [1 0 x, =y, || Ay| a=(1+m)cos6 tes 
yj, [0 1 y, x fla | b=(1+m)sing 
误差 方程 为 
Ax 
Va 1 0 x, y, lA Xp —X 
-| : g y -| : | a 
a 0 1 y x, la Yy = 
b 


A 
其 中 ， H 为 转换 前 坐标 ， H 为 转换 后 坐标 ， | | 为 平移 参数 ，(1+m) 为 缩放 参数 ， 
0 为 旋转 参数 。 理 论 上 ， 只 要 存在 两 组 转换 控制 点 ， 利 用 上 式 方程 ， 可 联 立 解 出 一 组 转换 
四 参数 ;存在 3 个 或 3 个 以 上 的 控制 点 时 ,可 利用 最 小 二 乘法 求解 四 参数 X = (Ax, Ay, O) ; 


然后 将 待 匹 配 图 像 转换 到 参考 图 像 。 

需要 注意 ， 在 视 场 平面 坐标 转换 时 ， 需 要 对 误 配 点 对 进行 筛选 ， 通 过 计算 控制 点 间 的 距 
离 。 在 实际 匹配 中 ， 如 果 控 制 点 之 间 的 距离 大 于 5 像 元 ， 即 2. 5” 时 ， 则 认为 该 点 为 误 配 点 
对 ， 将 其 剔除 ， 根 据 筛 选 后 的 控制 点 求解 的 RST 参数 ， 重 新 对 SDO/HMI 和 1M 太阳 望远镜 像 
进行 RST ISN, iHi RST 参数 ， 最 后 得 到 友 代 完成 后 的 RST 转换 参数 。 


3 实验 结果 与 分 析 


本 文 实验 中 ， 分 块 尺寸 设置 为 30 像 元 *30 像 元 ， 块 间隔 为 5 像 元 。 图 3 为 利用 上 述 方 
法 ， 将 IM 太阳 望远镜 高 分 辨 像 和 SDO/HMI 全 日 面 像 配 准 后 重合 区 域 的 对 比 图 ， 


a IM 太阳 望远镜 像 b SDO/HMI 对 应 的 区 域 
图 3 1M 太阳 望远镜 高 分 辨 像 和 SDO/HMI 全 日 面 像 配 准 后 重 闭 区 域 的 对 比 图 


Fig. 3 Comparison overlapping area of 1M 太阳 望远镜 high-resolution image and full-disk Solar image after 


registration 
从 图 3 可 以 看 到 , 配 准 后 图 像 的 对 应 区 域 除 清晰 度 明显 不 同 外 , 其 他 特征 已 经 完全 匹配 。 
为 了 更 好 地 观察 配 准 效果 ， 将 配 准 区 域 看 加 后 ， 对 边界 区 域 进行 放大 ， 如 图 4， 可 以 看 到 边 
界 细节 的 信息 匹配 度 良 好 ， 和 气孔、 米粒、 米粒 暗 径 均 较 好 的 衔接 。 
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图 4 配 准 区 域 的 边界 细节 


Fig.4 The edge details of the registration area 


根据 配 准 结果 和 SDO/HMI 的 标准 参数 ， 表 1 显示 了 推算 的 这 张 IM 太阳 望远镜 图 像 PO 
角 和 像 元 比例 尺 ， 以 及 视 场 中 心 在 SDO/HNI 全 日 面 像 中 的 位 置 。 


表 1 配 准 结果 


Tab. 1 Registration results 


PET 
POAR BIR som FAB ne, Beto 
配 准 结 CE) (arcsecond 像 元 ) ( 像 元 ) 标 匹配 点 l 数 (SDO 像 元 ) 
终 配 准 结 (ez) 
R 
-12. 95° 0. 04223 1524. 163 2686. 251 775 0. 31 像 元 


从 表 1 可 以 知道 平均 残 差 为 0. 31SD0/MI 像 元 ， 意 味 着 拟 合 结果 的 偏差 约 为 0. 15 ”。 
5 显示 了 这 些 残 差 的 空间 分 布 。 从 图 5 可 以 看 出 ， 残 差分 布 并 不 是 完全 随机 ， 显 示 了 1M 
太阳 望远镜 像 存 在 畸变 。 这 是 由 于 地 面 高 分 辨 观测 图 像 重 建 后 , 依然 还 有 一 些 低频 的 大 气 拌 


动 残余 未 能 消除 , 而 这 种 原因 限制 了 最 终 拟 合 结果 的 精度 。 同时 由 于 两 个 图 像 的 观测 时 间 并 


不 完全 同步 ， 有 6 s 的 间隔 ， 


因此 有 些 太阳 图 像 特征 已 经 发 生 了 改变 。 
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图 5 配 准 后 的 向 量 场 


Fig.5 The vector field after registration 


对 多 张 1M 太 阳 望 远 镜 和 新 太阳 望远镜 观测 的 高 分 辨 像 和 SDO/HMI 全 日 面 像 进行 了 匹配 ， 
结果 见 表 2。 
表 2 新 太阳 望远镜 和 1M 太阳 望远镜 的 Ti0 图 像 与 SD0/HMI 全 日 面 像 的 配 准 结果 
Tab. 2 Registration result of TiO image( 新 太阳 望远镜 and IM 太阳 望远镜 ) and SDO/HMI 
full-disk Solar image 


5 SDO/HMI 的 残 差 均 方 根 


序号 望远镜 高 分 辨 像 观测 时 间 ine 2 
观测 间隔 (SDO 像 元 ) 
1 新 太阳 望远镜 2014-08-05 22:29:59 ls 0. 33 
2 新 太阳 望远镜 2014-08-05 16:50:13 2s 0. 43 
3 新 太阳 望远镜 2014-08-01 16:53:16 ls 0. 45 
4 新 太阳 望远镜 2012-07-05 16:36:28 Ts 0. 42 
1M 太阳 望 远 
5 2013-07-15 07:29:09 6s 0.31 
镜 
1M 太阳 望 远 
6 2012-10-29 06:03:37 8s 0. 50 
镜 
1M 太阳 望 远 
7 " 2013-06-12 07:45:34 ll s 0. 26 
Bi 


从 上 述 结果 可 以 看 到 ， 最 大 的 匹配 残 差 约 0. 5 像 元 ， 相 当 于 0. 25” « 
实验 过 程 中 发 现 , 不 同 参数 的 选取 会 对 算法 产生 一 定 的 影响 ， 具 体 见 图 6- 图 8， 
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Residual / Pixel 


->- 30 degrees | 
-0- 35 degrees 
-*- 40 degrees 
-+ 45 degrees 7 


图 6 不 同 的 初始 参数 


Fig.6 The different initial parameters 
图 6 显示 了 初步 估计 的 参数 对 迭代 次 数 的 影响 。 图 6 中 , 三 角 型 (绿色 )、 空 心 圆 型 ( 蓝 
E) KFE AEN TFE ERE) 分 别 代表 估计 角度 为 30”、35”、40”、45” 时 的 


夫 代 过 程 。 可 以 看 出 ， 随 着 初 值 误 差 增 大 ， 为 了 修正 误差 使 结果 收敛 在 合理 的 范围 内 ， 算 法 


的 迭代 次 数 相 应 增加 。 当 初始 估计 和 角度 达到 45” 《此 时 的 误差 约 为 30” ) 时 ， 和 迭代 次 数 由 


6 次 增加 到 42 次 才能 收敛 。 


lterations 
a 8 BB 8 & è 
T —T 


a 38 
T 


-+- 5 Pixel 


-0- 10 Pixel | 


1 1 1 1 1 
20 25 30 35 40 
Angle error / Degree 


图 7 不 同 的 控制 点 筛选 间距 


Fig.7 The different filter distance between control points 
图 7 显示 了 筛选 控制 点 间距 离 的 设置 对 算法 收敛 速度 的 影响 ,图 7 中 , 十 字 型 (红色 )、 
空心 圆 型 (黑色 ) 分 别 代表 筛选 阔 值 为 5 Pixel 和 10 Pixel 的 迭代 过 程 。 可 以 看 出 ， 初 值 


误差 较 大 时 ， 高 阔 值 的 迭代 次 数 更 少 ， 意 味 着 算法 收敛 更 快 ， 反 之 低 阔 值 效果 更 理想 。 当 角 


度 误差 在 30” 附 近 时 ,不 同 闪 值 的 迭代 次 数 相 差 28 次 , 高 闷 值 更 优 。 角度 误差 在 10” 附 近 


IY, ERARA OK, REL HE 
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-*- Sigma=0.5 
-0- Sigma=1.0 
-+- Sigma=1.5 


Residual / Pixel 
o o 


图 8 不 同 的 平滑 参数 


Fig.8 The different filtering parameters 


图 8 显示 了 不 同 的 平滑 参数 对 迭代 次 数 的 影响 。 图 8 中 , 米 字 型 (绿色 )、 空 心 圆 型 ( 蓝 


色 )、 十 字 型 〈 红 色 ) 分 别 表示 平滑 参数 为 0.5、1、1. 5 时 的 迭代 过 程 。 可 以 看 出 ,平滑 参 
数 增 大 使 迭代 次 数 增长 ， 并且 拟 合 结果 有 一 定 的 降低 。 当 平滑 参数 为 0.5 时 ， 算 法 达到 收敛 
所 需 的 迭代 次 数 和 结果 较 优 于 其 它 情况 。 


4 结论 


本 文 提出 了 一 种 针对 太阳 高 分 辨 像 和 全 日 面 进 行 视 场 匹 配 的 算法 , 从 而 解决 了 不 同 空间 


分 辨 率 下 、 具 有 旋转 角 的 太阳 观测 图 像 匹 配 问 题 。 算 法 的 基本 思想 是 结合 区 域 相关 匹配 和 控 


制 点 最 小 二 乘 求 解 的 方法 。1M 太阳 望远镜 和 新 太阳 望远镜 高 分 辨 图 像 与 SD0/HMI 全 日 面 像 
的 视 场 匹 配 实验 都 证 明了 两 者 的 图 像 拟 合 结果 的 偏差 优 于 0. 25”。 通 过 匹配 可 以 精确 测量 


高 分 辨 像 的 PO 角 ， 像 元 比例 尺 等 。 


这 套 方法 要 计算 大 量 的 小 区 域 相关 ， 同 时 进行 迭代 计算 因此， 时间 开 销 较 大 ; 以 及 不 


同 参数 对 和 妈 代 次 数 的 影响 ,是 在 方法 应 用 时 需要 注意 的 。 除 此 之 外 ,小 太 寸 子 块 会 导致 匹配 
点 数 增 加 ， 计 算 量 增 大 ， 拟 合 精度 会 存在 一 定 的 下 降 ， 而 大 斥 二 虽然 对 精度 有 一 定 提 升 ， 但 
是 会 导致 匹配 区 域 变 小 , 从 而 会 影响 整 张 图 像 的 匹配 效果 。 而 子 块 间距 更 多 的 影响 计算 速度 。 


结合 实际 情况 ， 我 们 认为 子 块 尺寸 在 30 像 元 、 块 间距 5 像 元 附近 ， 算 法 的 稳定 性 较 合理 。 
在 应 用 中 ， 可 根据 实际 情况 对 参数 进行 调整 。 


此 外 , 将 高 分 辨 Ti0 像 通过 与 SD0/HMI 全 日 面 连续 谱 图 像 的 高 精度 匹配 , 也 意味 着 可 以 


转换 为 和 任意 标准 的 全 日 面 像 进行 匹配 ， 如 SDO/HMI 磁 图 、SD0/AIA 各 波段 观测 ， 甚 至 和 色 
球 全 日 面 像 进 行 匹 配 。 这 意味 着 可 以 在 0.25 ”的 偏差 下 ， 研 究 小 尺度 结构 在 各 个 波段 ， 以 


及 在 磁场 上 的 对 应 情况 。 


这 套 方案 可 以 用 于 任何 具有 相似 图 像 结构 的 太阳 图 像 匹 配 , 如 Ti0 和 Gband, Ha 偏 带 、 


10830; 也 可 以 用 于 Ha 线 心 图 像 和 GONG 全 日 面 像 的 配 准 。 因 此 具有 非常 多 的 应 用 前 景 。 
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High-accuracy Registration Method of Solar High-resolution 


Observation Images and Full-disk Solar Images 


Feng Tao!, Dai Wei! , Wang Rui*, Ji Kaifan!? 

1.Key Laboratory of Applications of Computer Technology of the Yunnan Province, Kunming University of Science and Technology, 
Kunming 650500, China; 
2.Yunnan Observatories, Chinese Academy of Sciences, Kunming 650011, China 
Abstract: In the research of the Solar observation, the registration between high-resolution 
observation images and full-disk Solar images is a greatly meaningful work. However, it is 
difficult to register with high accuracy due to the rotation, scaling and translation between them. 
This paper presents an image registration method combining local statistical information and 
control point matching: At first, between the high-resolution observation image and the full-disk 
Solar image is preliminary estimated about the orientation, scale, and positional parameters. Then 
the images are pretreated based on the estimated parameters. And then the pretreatment of the field 
of view is divided into a large number of overlapping local regions equally. Each corresponding 
local area is determined by the Correlation matching on the full-disk after preprocessing. Next, the 
sub-pixel offset between each pair of local regions is then measured, and the coordinate positions 
of the feature point are determined according to the sub-pixel offset between each pair of local 
regions, which is used as the feature control point in the point matching. Finally, the conversion 
equations of the affine transformation are established according to the control point, then use the 
least squares to solve the entire field of view of the transformation parameters. The image is 
reiterated according to the parameters to be solved, and the transformation parameters are obtained 
after the final iteration to register images. By registering the high-resolution observation images 
and the full-disk Solar SDO/HMI continuous spectrum, the EMMOM onm is within 
0.25 arc-seconds. 


Key words: Solar Image Registration; High-resolution Image; Full-disk Solar Images; 


